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Abstract 

We give a review and present a comprehensive calculations for the leptonic 
constant /s c of the low-lying pseudoscalar and vector states of 5 c -meson in 
the framework of static and QCD-motivated nonrelativistic potential models 
taking into account the one-loop and two-loop QCD corrections in the short 
distance coefficient that governs the leptonic constant of B c quarkonium sys- 
tem. Further, we use the scaling relation to predict the leptonic constant of 
the nS-states of the be system. Our results are compared with other models 
to gauge the reliability of the predictions and point out differences. 

I. INTRODUCTION 

The B c mesons provide a unique window into the heavy quark dynamics. Although they 
are intermediate to the charmonium and bottomonium systems the properties of B c mesons 
are a special case in quarkonium spectroscopy as they are the only quarkonia consisting of 
heavy quarks with different flavours. Because they carry flavour they cannot annihilate into 
gluons so are more stable and excited £? c -states lying below BD (and BD* or B*D) threshold 
can only undergo radiative or hadronic transitions to the ground state B c which then decays 
weakly. There are two sets of S-wave states, as many as two P-wave multiplets (the IP and 
some or all of the 2P) and one D-wave multiplet lying below BD threshold for emission of B 
and D mesons. As well, the F-wave multiplet is sufficiently close to threshold that they may 
also be relatively narrow due to angular momentum barrier suppression of the Zweig allowed 
strong decays. However, the spectrum and properties of these states have been calculated 
various times in the framework of heavy quarkonium theory [1]. 

The discovery of the B c meson by the Collider Detector at Fermilab (CDF) Collaboration 
[2] in the channel B c — > J/i^lv (I = e, //, r) with low-lying level of pseudoscalar mass M Bc = 
6.4^o!i3 GeV and lifetime tb c = 0.46to;}g ± 0.03 ps confirmed the theoretical predictions 
regarding various B c meson properties, spectroscopy, production and decay channels [1]. 
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The description of various aspects of B c meson physics [3] has received recently a lot of 
attention both from theoretical as well as experimental sides. It could offer a unique probe to 
check the perturbative QCD predictions more precisely and lead to a new information about 
the confinement scale inside hadrons. The higher frc-meson state cascades down through 
lower energy be state via hadronic or electromagnetic transitions to the pseudoscalar ground 
state B c meson. The theoretical study of the pure leptonic decays of the B c meson such 
as B c — > Ivi can be used to determine the leptonic decay constant fs c [4], as well as the 
fundamental parameters in the standard Model (SM). Hence, it is the only meson within 
SM which is composed of two nonrelativistic heavy quarks of different flavors: open charm 
and bottom quarks. Its spectroscopy production mechanisms and decays differ significantly 
from those of charmonium J/ip and upsilon T as well as hadrons with one heavy quark. 
Indeed, this meson is also a long lived system decaying through electroweak interactions. It 
stands among the families of cc and bb and thus could be used to study both quantitatively 
and conceptually existing effective low energy frameworks for the description of bound state 
heavy quark dynamics, like NRQCD [5], pNRQCD [6] and vNRQCD [7]. 

The calculation of the B c leptonic decay constant /# c can be carried out either using 
QCD-based methods, such as lattice QCD [8], QCD sum rules [9], or adopting some con- 
stituent quark models [1]. So far, lattice QCD has only been employed to calculate the 
B c purely leptonic width. As the QCD sum rules, the B c leptonic constant, as well as the 
matrix elements relevant for the semileptonic decays were computed. Moreover, the leptonic 
constant was also been calculated by using the nonrelativistic (NR) potential model [1]. 

Ikhdair et al. [10] applied the NR form of the statistical model to calculate the spec- 
troscopy, decay constant and some other properties of the heavy mesons, including the be 
system, using a class of three static quarkonium potentials. Davies et al. [8] predicted the 
leptonic constant of the lowest state of the B c system with the recent lattice calculations. 
Eichten and Quigg [1] gave a more comprehensive account of the decays of the B c system 
that was based on the QCD- motivated potential of Buchmiiller and Tye [11]. Gershtein et 
al. [1] also published a detailed account of the energies and decays of the B c system using 
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a QCD sum-rule calculations [9]. Fulcher [2] also calculated the leptonic constant of the B c 
quarkonium system using the treatment of the spin-dependent potentials to the full radiative 
one-loop level and thus included effects of the running coupling constant in these potentials. 
Furthermore, he also used the renormalization scheme developed by Gupta and Radford [12]. 
Ebert et at. [1] have also calculated the leptonic constant of the pseudoscalar and vector B c 
quarkonium system using a relativistic model and then compared their calculation with the 
NR model. Kiselev et al. [13-15] calculated the leptonic constant of the pseudoscalar meson 
B c in the framework of QCD-motivated potential models taking into account the one-loop 
[13], the two-loop [14] correction matching and scaling relation (SR) [15]. Capstick and 
Godfrey [16] predicted the result of f Bc = 410 MeV for the decay constant of the B c meson 
using the Mock-Meson approach or other relativistic quark models. Ikhdair and Sever [17] 
have calculated the decay constant in the nonrelativistic and semi-relativistic quark model 
using the shifted large-N expansion technique. Finally, Godfrey [18] has calculated the spec- 
troscopy of B c meson in the relativixzed quark model. Motyka and Zalewski [18] proposed a 
nonrelativistic Hamiltonian with plausible spin dependent corrections for the quarkonia be- 
low their respective strong decay threshold. They found the decay constants for the ground 
state and for the first excited S-state of be quarkonium to be f Bc = 435 MeV and fs c =315 
MeV, respectively. 

We have made this study mostly a fully treatment for the quarkonium potential models 
used in the literature [1,17] in calculating the B c leptonic constant based on the NR potential 
model, (cf. e.g. [17] and references contained there), and applied with a great success to fit 
the entire heavy quarkonum spectroscopy [17,19]. We insist upon strict flavor-independence 
of its parameters. We also use the potential models to give a simultaneous account of the 
properties of the cc, bb and be quarkonium systems. 

The contents of this paper is as follows: In section II, we present briefly the solution of 
the Schrodinger equation for the be quarkonium mass spectrum using the shifted large-N- 
expansion technique (SLNET). In section III, we introduce the necessary expressions for the 
one-loop and two-loop corrections to the leptonic constant available to moment. Further, 
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in section IV we present the phenomenological and the QCD-motivated potentials used in 
the present work. Finally, section V contains our calculations and remarks made for the 
leptonic constant in our approach and scaling relation. 



II. WAVE EQUATION 



In this section we consider the N-dimensional space Schrodinger equation for any spherically 
symmetric potential V(r). If ip(r) denotes the Schrodinger's wave function, a separation of 
variables ip(r) = Y ljfn (9, (ft)u(r) /r^' 1 ^ 2 gives the following radial equation (in units k — c — 
1) [17,19,20] 

f 1 d 2 [fc - (1 - a)]\k- (3 - a)] 



+ 



+ V(r)\u(r)=E n , lU (r), 



(1) 



[ 4/i dr 2 16/zr 2 , 

where ji = (rn^rrig^j /(m^ + m q2 ) is the reduced mass for the two interacting particles. 
Here E n j denotes the Schrodinger binding energy and k = N + 21 — a, with a representing a 
proper shift to be calculated later on and / is the angular quantum number. We follow the 
shifted 1/N or 1/k expansion method [17] by defining 



Q V dr o J rn\ 



(2) 



and also the energy eigenvalue expansion [17,19] 



m=0 



Q 



-E 

i- 1 mi 



(3) 



—1/2 

where x = k (r/r — 1) with r is an arbitrary point where the Taylor expansion is being 
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performed about and Q is a scale to be set equal to k at the end of our calculations. 
Inserting equations (2) and (3) into equation (1) gives 



Id 2 1/k (2 - a) (l-a)(3-a) 



(1 



Andy 2 Ajj, \ 4 



4k 



E (-i) r 



m) 



m=0 



r 2 oo 

+- E 



'd m V(r ) 
dr^ 



NTp(m-2)/2 



ml 



Xn r 



(4) 
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where the final analytic expression for the 1/k expansion of the energy eigenvalues appro- 
priate to the Schrodinger particle is [17] 



r 2 00 n \ 
f - !l V k ( m) F 

r ~ Q 



(5) 



m=0 



where n r is the radial quantum number. Hence, we formulate the SLNET (expansion as 1/k) 
for the nonrelativistic motion of spinless particle bound in spherically symmetric potential 
V(r). The resulting eigenvalue of the N-dimensional Schrodinger equation is written as [17] 



in r — k 



1 + r 2 V(r ) 



1 



I6/1 Q 
1 — a) (3 — a 



+ 



t 1 , (2 -a) 
(n r + ± '- 

K r T 8/1 



16// 



+ a 



(i) 



a 



(2) 



(6) 



where and are two useful expressions given by Imbo et al [20]. Comparing equation 
(5) with (6) yields 

Q 



E = V(r ) + 



16 fir 2 



2 • 



(7) 



Ei — Q 



1\ (2 -a) 

n r H — uj 

2) 8/i 



(8) 



E ~ Q 



:i-«)(3-a) + a w 



and 



16 fi 



Q 



(9) 



E 3 = 4a< 2 >. 



r 2' 



(10) 



The quantity r is chosen so as to minimize the leading term, E [17,19], that is, 

dE d 2 E 

■ and — tt > 0. 



dr, 



Therefore, r satisfies the relation 



drl 



(11) 



Q = &Hr'V'{ro), 



(12) 
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and to solve for the shifting parameter a, the next contribution to the energy eigenvalue E\ is 
chosen to vanish [20] so that the second- and third-order corrections are very small compared 
with the leading term contribution. The energy states are calculated by considering the 
leading term E and the second-order and third-order corrections, it implies the shifting 
parameter 

r V"(r y 1/2 



a = 2- (2n r + 1) 



3 + 



V'(r ) 



(13) 



Therefore, the Schrodinger binding energy to the third order is 



B „,, = nr . o) + !^ + i 

Z ' o 



'1 — a) (3 — a) m a 



(2) 



(14) 



16/x k \k' 

Once the problem is collapsed to its actual dimension (N = 3), it simply rests the task 
of relating the coefficients of our equation to the one-dimensional anharmonic oscillator in 
order to read the energy spectrum. For the N = 3 physical space, the equation (1) restores 
its three-dimensional form. Thus, with the choice k = \[Q which rescales the potential, we 
derive an analytic expression that satisfies r in equations (12)-(13) as 



1 + 21 + (2n r + 1) 



3 + 



V'(r ) 



8KV>o) 



1/2 



(15) 



where n r = 0, 1, 2, • • • stands for the radial quantum number and / = 0, 1, 2, • • • stands 
for the angular quantum number. Once r is being determined through equation (15), the 
Schrodinger binding energy of the q x q2 system in equation (14) becomes relatively simple and 
straightforward. We finally obtain the total Schrodinger mass binding energy for spinless 
particles as 



^(<?i<?2) = + m q2 + 2E n> i. 



(16) 



As stated before in [17,19,20], for any fixed n the computed energies become more convergent 
as / increases. This is expected since the expansion parameter 1/N or 1/k becomes smaller as 
/ becomes larger since the parameter k is proportional to n and appears in the denominator 
in higher-order correction. 



On the other hand, the spin dependent correction to the nonrelativistic Hamiltonian, 
which is responsible for the hyperfine splitting of the lS-state mass level is generally used 
in the form (cf. e.g. Refs. [l^^l]) 1 

Like most authors (cf. e.g. [1,17,21]), we determine the coupling constant a s (m L ) from the 
well measured experimental hyperfine splitting of the lS(cc) state value 

A£ H f = M J/% b - M Vc = 117 ± 2 MeV. (18) 

We use the value of the coupling constant to reproduce the spin-averaged data (SAD) or 
center-of-gravity (c.o.g.) of the lowest charmonium state value, that is, M^(1S). In order to 
apply this formula one needs the value of the wave function at the origin, this is obtained 
by solving the Schrodinger equation with the nonrelativistic Hamiltonian and the coupling 
constant a s (/i). The radial wave function at the origin [21] is determined by 

|rf = 2„(^). (19) 

where |i?is(0)| 2 = 4n \ipis(0)\ 2 . Hence, the total mass of the low-lying pseudoscalar £? c -state 
is given by the expression 

M Bc ((T) =m c + m b + 2E lfi - 3AE UF /4, (20) 

and also for the vector £?*-state 

Mb c .(1 _ ) =m c + m b + 2E lfi + AE UF /A. (21) 

Finally, the square-mass difference can be simply found by using the expression 

AM 2 = M|»(l~) - M§ (0~) = 2AE H f [m c + m b + 2E lfi - AE H f/4] . (22) 



l To the moment, the only measured splitting of nS- levels is that of rj c and J/ip, which allows us 
to evaluate the so-called SAD using M^(1S) = (3M J/V , + Mr] c )/4 and also M(nS) = M v (nS) - 
(Mj^-Mri c )/4n [13,21]. 

8 



III. LEPTONIC DECAY CONSTANT OF THE 5 c -MESON 



The study of the heavy quarkonium system has played a vital role in the development of the 
QCD. Some of the earliest applications of perturbative QCD were calculations of the decay 
rates of charmonium [22]. These calculations were based on the assumption that, in the NR 
limit, the decay rate factors into a short-distance (SD) perturbative part associated with 
the annihilation of the heavy quark and antiquark and a long-distance (LD) part associated 
with the quarkonium wavefunction. Calculations of the annihilation decay rates of heavy 
quarkonium have recently been placed on a solid theoretical foundation by Bodwin et al. 
[23]. Their approach is based on NRQCD, and effective field theory that is equivalent to 
QCD to any given order in the relative velocity v of the heavy quark and antiquark [24]. 
Using NRQCD [5] to seperate the SD and LD effects, Bodwin et al. derived a general 
factorization formula for the inclusive annihilation decay rates of heavy quarkonium. The 
SD factors in the factorization formula can be calculated using pQCD, and the LD factors 
are defined rigorously in terms of the matrix elements of NRQCD that can be estimated 
using lattice calculations. It applies equally well to S-wave, P-wave, and higher orbital- 
angular-momentum states, and it can be used to incorporate relativistic corrections to the 
decay rates. 

In the NRQCD [25] approximation for the heavy quarks, the calculation of the leptonic 
decay constant for the heavy quarkonium with the two-loop accuracy requires the matching 
of NRQCD currents with corresponding full-QCD axial-vector currents [14] 

J x = -xl*p c v x and J A = &7 A 7 5 c, (23) 

NRQCD A6rc QCD ' /0 ' V ' 

where b and c are the relativistic bottom and charm fields, respectively, xl an d 4> c are the 
NR spinors of anti-bottom and charm and v A is the four- velocity of heavy quarkonium. The 
NRQCD [25] lagrangian describing the £> c -meson bound state dynamics is 

£nrqcd = Aight + ^ (iD + D 2 /(2m c )) ^ + xl (iA> - D 2 /(2m 6 )) Xb + • • • , (24) 
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where Cu g u is the relativistic lagrangian for gluons and light quarks. The two-component 
spinor field ip c annihilates charm quarks, while Xb creates bottom anti-quarks. The relative 
velocity v of heavy quarks inside the £? c -meson provides a small parameter that can be used 
as a nonperturbative expansion parameter. To express the decay constant fs c in terms of 
NRQCD matix elements we express J A in terms of NRQCD fields ip c and Xb- The A = 



QCD 



current-component contributes to the matrix element and consequently the J x 
following operator expansion 



has the 

QCD 



(0\b 1 x l5 c\B c (P))= i f B P\ (25) 

where \B C (P)) is the state of the 5 c -meson with four-momentum P. Only the A = compo- 
nent contributes to the matrix element (25) in the rest frame of the 5 c -meson. It has the 
standard covariant normalization 

(2^)8 / r Bc (P^B c (p)d 3 P = 2£5 (3 V - p), (26) 

and its phase has been chosen so that f Bc is real and positive. In NRQCD, the matching 
yields 

6 7 75 c = K xtiJ c + K 2 (B Xb y.-Dij c + • • • , (27) 

where K = K (m c ,m b ) and K 2 = K 2 (m c ,mb) are Wilson SD coefficients. They can be 
determined by matching perturbative calculations of the matrix element (0| b^'y^c \ B C ) that 
is mainly resulting from the operator xl^c i n 



5 7 °7 5 c \B C ) ^ = K (0| xhc \B C ) mnnn + K 2 (0| (D X6 )t.D^ c \B 



QCD 



+ (28) 

NRQCD 



NRQCD 

where the matrix element on the left side of (28) is taken between the vacuum and the state 
\B C ) . Hence, equation (28) can be estimated as 

<0|x^ c |£ c > 2 ^^|^is(0)| 2 . (29) 

71 

Onishchenko and Veretin [25] calculated the matrix elements on both sides of equation (28) 
up to a 2 s order. In one-loop calculation, the SD-coefncients are 
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K = 1 and K 2 = 



8/x 2 ' 



(30) 



with /j, denned after equation (1). Further, Braaten and Fleming in their work [26] calculated 
the perturbation correction to K up to order- a s (one-loop correction) as 



Kr 



with ci being calculated in Ref. [26] as 



Cr 



7Y 



(31) 



ci = - 



m b -m c m b 
2 in — 



777.6 + m c m c-i 

Finally, the leptonic decay constant for the one-loop calculations is 



(32) 



(l-loop) 
B c 



1 - 



a s (fi) 



7Y 



m b - m c m b 
2 in — 

777;, + fn c m c 



fNR 

JB C ) 



(33) 



where the NR leptonic constant [27] is 



NR 



l^is(0)| 



(34) 



and fj, is any scale of order 777 fe or 777 c of the running coupling constant. On the other hand, 
the calculations of two-loop correction in the case of vector current and equal quark masses 
was done in [28] . Further, Onishchenko and Veretin [25] extended the work of Ref. [28] into 
the the non-equal mass case. They found an expression for the two-loop QCD corrections 
to _B c -meson leptonic constant given by 



K (a s , M/n) = 1 + Cl (M//i)^S + C2 (M//i) 

7T \ 71 



+ 



(35) 



where C2(M/fj) is the two-loop matching coefficient and with are explicitly given in Eq. 
(32) and (Ref. [25]; see equations (16)-(20) therein), respectively. In the case of £? c -meson 
and pole quark masses (m b = 4.8 GeV, 777 c = 1.65 GeV), they found 



(2-loop) 



1-1.48 f^^- 24.24 



71 



71 



fNR 
JB C ■ 



(36) 



Therefore, the two-loop corrections are large and constitute nearly 100% of one-loop correc- 
tion as stated in Ref. [25]. 
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IV. SOME POTENTIAL MODELS 

The potential models are suitable for the phenomenological studies, because they can re- 
produce the model formulae or numbers for the quantities, used as input values (the level 
masses, for example). Hence the potential models can be considered as phenomenological 
meaningful fittings of some experimental values, but they can not restore a true potential, 
that does not exist due to the nonperturbative effects [15]. 

A. Static potentials 

It is easy to see that most phenomenological static potentials used in the Scrodinger 
equation in [17] may be gathered up in a general form [17,29] 

V(r) = -ar- a + br + V 0<a,(3<l, a,6>0, (37) 

where V may be of either sign. Here we assume that the effective potential consists of 
two terms, one of which, Vy{r) = —ar~ a , transforms like a time-component of a Lorentz 4- 
vector and the other, Vs(r) = br 13 , like a Lorentz scalar. These static quarkonium potentials 
are monotone nondecreasing, and concave functions which satisfy the condition 

V\r) > and V"{r) < 0. (38) 

At least ten potentials of this general form, but with various values of the param- 
eters, have been proposed in the literature (cf. [17] and references contained there). 
The generality (37) comprises the following five types of potentials used in the litera- 
ture: Some of these potentials have a = (3 as in the Cornell (a — f3 — 1); Song- Lin 
(a — P — \) and Turin (a — (5 — |) potentials with same sets of fitting parame- 
ters used in our previous works [17]. Further, Song [30] has also used a potential with 
a — (5 — |. 

On the other hand, potentials with a ^ (3 have also been popular as 
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1. Martin potential 



The phenomenological power-law potential proposed by Martin (cf. e.g. [15,17]) has 
a — 0, P — 0.1 of the form 

V M (r) = b M (A M r)°- 1 + c M , (39) 
is labeled as Martin's potential [15] with a given set of adjustable parameters 

[6m, cm, Am] = [6.898 GeV 11 , -8.093 GeV, 1 GeV] , (40) 

and quark masses 

[m c , m b ] = [1.800 GeV, 5.174 GeV] . (41) 
(potential units are also in GeV). 

2. Logarithmic -potential 

A Martin's power-law potential turns into the logarithmic potential of Quigg and Rosner 
[15,17] corresponds to a = [3 — > and it takes the general form 

VL(r) = 6 L ln(A L r) + c L , (42) 

with an adjustable set of parameters 

[b u c L , A L ] = [0.733 GeV, -0.6631 GeV, 1 GeV] , (43) 

and quark masses 

[m c , m b ] = [1.500 GeV, 4.905 GeV] . (44) 

The potential forms in (39), and (42) were also used by Kiselev in [15]. Further, they were 
also been used for ip and T data probing 0.1 fm < r < 1 fm region [13]. 

Further, Motyka and Zalewski [29] used a nonnrelativistic potential with a = 1 and 
(3 — \ for the bb quarkonium and then applied it to the cc and be quarkonium systems. 
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Grant, Rosner, and Rynes [31] have suggested a = 0.045, (3 — 0. Heikkila, Tornquist and 
Ono [32] tried a — 1, (3 — |. Some very successful potentials known from the literature are 
not of this type. Examples are the Indiana potential [33] and the Richardson potential [34]. 

B. QCD-motivated potentials 

We use two types of the QCD-motivated potentials: The Igi-Ono (10) and an improved 
Chen-Kuang (CK) potential models. The details of these potentials can be tracesd in our 
previous works in [17]. 

V. CALCULATIONS AND CONCLUSION 

Based on our previous works [17], we determine the position of the charmonium center- 
of-gravity (c.o.g.) M^(1S) mass spectrum and its hyperfine splittings by fixing the coupling 
constant a s {m c ) in (17) for each central potential. Further, we calculate the corresponding 
low-lying (c.o.g.) Mx(lS) and consequently the low-lying Mb c (1S). 2 In Table I, we estimate 
the radial wave function of the low-lying state of the be system, so that 

|i?is(0)| = 1.18 - 1.24 GeV 3/2 , (45) 

for the set of the central potentials given in section IVA. Further, in Table II we present our 
results for the NR leptonic constant f^f = 466 ± 19 MeV and f^f = AQA ± 19 MeV as an 
estimation of the potential models without the matching. Our results are compared with 
those of Gershtein, Likhoded and Sabespi [35], who used Martin's potential and with those 
of Jones and Woloshyn [36]. Moreover, the one- loop correction /b* -1o ° p ^ and the two- loop 
correction / , £ c ~ loop ' ) are also given in Table III. Therefore, in the view of our results, our 
prediction for the one-loop calculation is 



2 Kiselev et al. [15] have taken into account that AM T (1S) = ^U|AM V ,(1S) with a s (T) / a s (ip) ~ 
3/4. On the other hand, Motyka and Zalewiski [29] also found ^ 11/18. 
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jO-loop) = 40 g ± 16 Me y ? ( 4 g) 

and also for two-loop calculation 

yg-ioop) = 315 ± 50 MeV. (47) 

Our estimation for /|J R is fairly in good agreement with the estimates in the framework 
of the lattice QCD result [8], = 440 ± 20 MeV, QCD sum rules [9], potential models 
[2] and also the scaling relation [15]. It indicates that the one-loop matching in Ref. [26] 
provides the magnitude of correction of nearly 12%. However, the recent calculation in the 
heavy quark potential in the static limit of QCD [13] with the one-loop matching [14] is 

jO-loop) = 40 q ± 15 Me y ; ( 4 g) 

Therefore, in contrast to the discussion given in [14], we see that the difference is not crucially 
large in our estimation to one-loop value in the _B c -meson. 
Our final result for the two-loop calculations is 

fV°° P) = 315±|g MeV, (49) 

the larger error value in (49) is due to the strongest running coupling constant in Cornell 
potential. 

Slightly different additive constants is permitted in this sector for a charmed mesons 
to bring up data to its (c.o.g.) value. However, with no additive constant to the Cornell 
potential [37], we notice that the smaller mass values for the composing quarks of the meson 
leads to a rise in the values of the potential parameters which in turn produces a notable 
lower value for the leptonic constant as seen in Tables II and III. 

On the other hand, for the QCD motivated two types 10 potential, our calculations 
for the ground state radial wave functions are presented in Table IV. In Table V, the NR 
leptonic constant of the pseudoscalar B c state are in the range 354 — 426 MeV, and for the 
B* are 353 - 424 MeV, for the type I. Further, they are found 353 - 457 and 351 - 454, for 
the type II. For instance, we may choose /|J R = 391 MeV with = 200 MeV, for the 
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type I, and /£ c R = 396 MeV with = 300 MeV, for the type II. These results surpass 
the matching procedure as well as methods followed in other references, cf. Refs. [8,13,14]. 
Further, for the CK potential, we present the decay constant in Table VI. 
The scaling relation (SR) for the S-wave heavy quarkonia has the form [15] 

f 2 n _ ( M n (bc) \ 2 f m c + m b \_d 
M n (bc) {M^bc)) \ 4// ) n' 1 ' 

where m c and m& are the masses of heavy quarks composing the _B c -meson, \i is the reduced 
mass of quarks, and d is a constant independent of both the quark flavors and the level 
number n. The value of d is determined by the splitting between the 2S and IS levels or 
the average kinetic energy of heavy quarks, which is independent of the quark flavors and 
n with the accuracy accepted. The accuracy depends on the heavy quark masses and it is 
discussed in [15] in detail. The parameter value in (50), d ~ 55 MeV, can be extracted from 
the experimentally known leptonic constants of tp and T [15]. Hence, in the view of Table 
VII, the SR gives for the lS-level 

f ( l R) * 444+^3 MeV, (51) 

for all static potentials used. Kiselev estimated f Bc = 400 ± 45 MeV [14] and = 
385 ± 25 MeV [15], Narison found /g R) = 400 ± 25 MeV [38]. 

Overmore, we give the leptonic constants for the excited nS-levels of the be quarkonium 
system in Table VII. We remark that the calculated value of f^jls) = ±15 MeV is in 
good agreement with fs^zs) = 280 ± 50 MeV being calculate in [13] for the 2S-level in the 
be system and it is also consistent with the scaling relation [15] . 

In Figure 1, we plot the calculated values of B c leptonic constants using (50) for different 
potential models together with the calculated values of the excited nS-states using [15] 

fn = -Wl. (52) 

Hence, our findings obtained from the potential model and SR are in good agreement to 
several MeV as shown in Figure 1. 
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Finally, we have also noted that the leptonic constant is practically independent of the 
total spin of quarks, so that 

fv,n ^ fp,n = fn, (53) 

where M Bc ~ m b + m c . Thus, one can conclude that for the heavy quarkonium, the QCD 
sum rule approximation gives the identity of leptonic constant values for the pseudoscalar 
and vector states. 

In this work, we have successfuly applied the SLNET using a class of static and QCD- 
motivated potentials to calculate numerically the leptonic constant of the pseudoscalar and 
vector S c -meson. Once the experimental leptonic constant of the B c meson becomes clear, 
one can sharpen the analysis. 
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TABLES 



TABLE I. The characteristics of the wave function | ^is (0) | 2 and the radial wave function 
\R 1S (0)\ = 47r|V>is(0)| both at the origin (in GeV 3 ) obtained from the Schrodinger equation for 
various potentials. 



Level 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


CornelP 


l^is(0)| 2 


0.112 


0.123 


0.111 


0.119 


0.102 


0.065 


\Ris(o)\ 2 


1.413 


1.54 


1.397 


1.495 


1.28 


0.814 



a Here we cite Ref. [37] for the fitted set of parameters. 



TABLE II. Pseudoscalar and vector decay constants (/p = /b c , /y = fs*) of the B c meson, 
calculated in the different potential models (the accuracy is 5%), in MeV. 



Constants 


CornelP 


Song-Lin 


Turin 


Martin 


Logarithmic 


Cornell b 


[35] 


[36] 


fNR 
JB C 


464.5 


485.1 


462.0 


478.1 


441.7 


351.5 


460±60 


420±13 


fNR 

Jbi 


461.5 


482.2 


459.2 


475.3 


439.3 


349.7 


460±60 





a Here V / 0. 

b Here V = 0, see Ref. [37]. 



TABLE III. One-loop and two-loop corrections to pseudoscalar and vector decay constants of 
the low-lying be system, calculated in the different potential models (the accuracy is 4%), in MeV. 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


Cornell 


/•(l— loop) 
J B c 


393.6 a 


424.4 


399.6 


421.2 


399.3 


311.9 


j,(2— loop) 
J B c 


264.1 b 


333.0 


296.6 


339.2 


340.9 


238.1 


r(l — loop) 

Jb* 


391.0 


421.9 


397.1 


418.7 


397.1 


310.3 


r-(2 — loop) 

Jbi 


262.4 


331.0 


294.8 


337.2 


339.0 


236.9 



a To the first correction for all potentials, Kq = 0.85 — 0.90. 
b To the second correction for all potentials, Kq = 0.57 — 0.77. 
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TABLE IV. The characteristics of the radial wave function at the origin |i?is(0)| 2 (in GeV 3 ) 
obtained from the Schrodinger equation for the Igi-Ono potential model. 



Lcvcl a 


1UU 


200 


300 


400 


500 


Type I b 












Uju Q (V\\\ 2 


0.066 


0.08 


0.092 


0.095 


0.089 


\Rism 2 


0.826 


1.005 


1.156 


1.19 


1.114 


Type IF 












I^15(0)| 2 


0.065 


0.071 


0.082 


0.096 


0.109 


\Rism 2 


0.819 


0.891 


1.03 


1.204 


1.366 



a For the fitted parameters, see Ref. [17]. 

b We shifted the cc spectra by c = -22 to -31 MeV. 

c We shifted the cc spectra by c = —15 to —26 MeV. 
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TABLE V. Pseudoscalar and vector decay constants of the B c meson, calculated in the Igi-Ono 
potential model, in MeV. 









1V MS 






fill Qirl"l "("AT^ 


i nn 

1UU 




ouu 


inn 


^nn 

ouu 


lype i .a s — 


n 1 qq c 




U.ZoO 


n 9^ 
u.zo 


n 9R9 


Ib R 


354.1 


391.1 


420.0 


426.0 


411.7 


fNR 
JB* 


352.6 


389.2 


417.7 


423.6 


409.4 


Type ll d :a s = 


0.199 


0.227 


0.230 


0.241 


0.251 


fNR 
J B c 


352.7 


368.2 


396.1 


428.6 


456.9 


fNR 

Jb* c 


351.2 


366.4 


394.1 


426.4 


454.4 



a For the fitted parameters we cite Ref. [17]. 
b We shifted the cc spectrum by c = —22 to —31 MeV . 
c Input value AEhf= 117 MeV; determines a s (m c ). 
d We shifted the cc spectra by c = —15 to —26 MeV. 



TABLE VI. The radial wave function, pseudoscalar and vector decay constants of the B c meson, 
calculated in the improved Chen-Kuang potential model, in MeV. 



Quantity 3, 


100 


180 


300 


350 


375 


400 


450 


500 


\Ris(o)\ 2 


1.017 


1.017 


1.017 


1.017 


1.017 


1.015 


0.957 


0.829 


fNR 
JB C 


393.5 


393.5 


393.5 


393.4 


392.5 


389.9 


381.8 


353.8 


fNR 

JB' c 


391.4 


391.4 


391.4 


391.3 


390.4 


387.8 


379.8 


352.3 



a Input value AEhf= 117 MeV; determines a s (m c ) = 0.27. 
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TABLE VII. The nS-levels leptonic constant of the be system, calculated in the different po- 
tential models (the accuracy is 3 — 7%), in MeV, using the SR. 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


/is 


449.6 


450.4 


448.0 


448.8 


420.9 


/2s 


305.8 


305.0 


303.3 


303.5 


284.7 


/3s 


243.0 


243.2 


241.3 


241.8 


227.2 


fas 


206.0 


207.1 


204.9 


205.9 


193.8 
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FIGURES 

FIG. 1. The nS-levels leptonic constant of the B c system calculated in the different potential 
models using SR. 
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